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Melanins are a group of complex pigments of biological origin, widely 
spread in all species from fungi to man. Among diverse types of melanins, 
the human melanins, eumelanins, are brown or black nitrogen-containing 
'b) . pigments, mostly known for their photoprotective properties in human skin. 

We have undertaken theoretical studies aimed to understand absorption 
spectra of eumelanins and their chemical precursors. The structure of the 
' biopigment is poorly defined, although it is believed to be composed of 

I cross- linked heteropolymers based on indolequinones. As a basic model 

. of the eumelanin structure, we have chosen pentamers containing hydro- 

' quinones (HQ) and/or 5,6-indolequinones (IQ) and/or semiquinones (SQ) 

often listed as structural melanin monomers. The eumelanin oligomers 
^Lj' have been constructed as random compositions of basic monomers and 

I opitimized for the energy of bonding. Absorption spectra of model assem- 

0^. blies have been calculated within the semiempirical intermediate neglect 

of differential overlap (INDO) approximation. Model spectrum of eume- 
lanin has been further obtained by sum of independent spectra of singular 
' polymers. By comparison with experimental data it is shown that the 

H I INDO/CI method manages to reproduce well characteristic properties of 

experimental spectrum of synthetic eumelanins. 

PACS numbers: 87.15.-v; 31.10.+z; 31.15.Gy; 36.20-r 
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1. Introduction 

Melanins are ubiquitous natural pigments formed by oxidation and poly- 
merization of catechols Black or dark brown nitrogen-containing 
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melanins are classified as eumelanins distinctly from pheomelanins - whose 
lighter coloring originates from intervention of cysteine in the process of 
melanogenesis, and allomelanins - which originate from nitrogen-free pre- 
cursors and are typical of plants and microorganisms. The melanogens - 
uncoloured precursors of melanins of the animal kingdom - are diphenols 
derived from phenylalanine and from tyrosine. The melanin pigments are 
amorphous substances of remarkable ability to absorb almost indiscrimi- 
nately near-infrared, visible and ultraviolet radiation [^|^. Exploiting these 
properties of melanins has proven to be a complex task, since despite sig- 
nificant experimental effort on a variety of natural and synthetic melanins, 
the molecular structure of the pigment and its organization remains still un- 
known. Poorly understood are also functional properties of natural melanin 
pigments which call for attention because of their both - photoprotective 
and photosensitizing properties. 

In this paper we will focus on hypothesis which concludes on spectroscopic 
properties of melanins assuming their polymeric amorphous structure. Our 
studies concern eumelanins which are known as complex broadband ab- 
sorbers, frequently subject to excitation in its native environments. Eu- 
melanin is mainly derived from 5,6-dihydroxyindole (DHI; for simplicity 
named also hydroquinone, HQ) and 5,6-dihydroxyindole-2-carboxylic acid 
(DHICA) j4|. The detailed structural properties of this polymer are still 
under study because of unclear covalent and ionic configurations of HQ 
and DHICA monomers within the natural eumelanin unit. It is possible, 
however, to produce synthetically pure DHI or DHICA melanin. In this 
respect, theoretical studies presented in this paper relate to synthetic eu- 
melanins that are free or almost free from the DHICA contribution in the 
structure. 

Absorption bands of eumelanins have been recently shown to depend 

on different levels of aggregation reported in scanning electron microscopy 
studies. The proposed aggregation structure is generally viewed as a molec- 
ular system containing 3-4 oligomers, often referred to in literature [HI C] 
as the fundamental aggregate. The hypothetical pigment structure is then 
believed to be assembled from these vr stacked, crosslinked units. 

In an earlier attempt to explain electronic (and optical) properties of melanins 
(c/. Figure 1), a band model [HUH] viewing melanins as amorphous semicon- 
ductors has been proposed. In fact, melanins can be considered as mixtures 
of more or less similar polymers, apparently made up of different structural 
units linked through heterogenous non-hydrozable bonds. The latter are the 
result of the polymerization of indole-type rings in both, the quinone and 
hydroquinone oxidation state, randomly crosslinked and piled. Theoretical 
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contributions on a speculative model 111 1 pictured eumelanins as a lin- 
ear chain semiconducting polymer. Extrapolating the bonding character of 
the lowest unoccupied molecular orbital (LUMO) of one particular dimer 
of 5, 6-indolequinone to the lowest conduction band of the infinite polymer, 
Pullman and Pullman have pointed out the tendency of such a melanin 
model to be an electron acceptor that could explain trapping of free radicals 
characteristic for the natural and synthetic pigments. 

More recent studies along similar lines have been performed by Galvao and 
Caldas |12l 1131 114j . Using the Hiickel vr-electron approximation and the 
parametrization of Pullman the authors studied the electronic struc- 
ture of a family of ideal ordered polymers arising from the indolequinone in 
different redox states. It has been shown that the direction of polymeriza- 
tion begins to emerge as the length of the polymer increases. The authors 
pointed out also that the redox state of the melanin units played a role 
in its band structure, as e.g. the polymer built from 5,6-dihydroxyindole 
units exhibited larger gaps and narrower bands, whereas finite chains of 
semiquinone units exhibited bonding character of their LUMO. 

Some limitation of this approach is the assumption of planarity of melanin 
polymers, which in most cases does not need to be the case. In fact, struc- 
tural modelling of eumelanins based on comparing calculated data of the 
reduced structure factor and radial distribution function for limited random 
network models with the experimental X-ray scattering data |15j suggest 
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Fig. 1. Experimental steady state absorption spectrum of synthetic DOPA melanin 
recorded using a UV-Vis spectrophotometer with an accuracy of 2nm. 
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that planar models consisting of undistorted long chains are not appropri- 
ate structures for amorphous melanins. 

We have thus considered the role of conformational variations of eumelanins' 
oligomers in modifying and modulation of their light- absorption properties. 
According to the existing evidence ^[^ElElEl that the 5,6-indolequinone 
(IQ) or the reduced forms, semiquinone (SQ) and hydroquinone (HQ) com- 
pose the major part of the active pigment, these three basic units (c/. Figure 
^) in their neutral forms have been used to build up model polymeric struc- 
tures. 

When focusing on oligomeric models of eumelanins, we have assumed in 
consistency with the wide-angle X-ray diffraction measurements of synthetic 
eumelanin samples and STM images of monolayers of synthetic eumelanins 
on graphite 0^1) that a fundamental building block of the polymer is a 
structure consisting of several (4-10) 5,6-indolequinone residues composed in 
TT-system molecular stacks. Therefore, as a starting model unit for a higher 
order polymeric structure of eumelanin a pentamer containing HQ, IQ or 
SQ molecules has been chosen and subsequently used to analyze randomly 
created structures of various conformations. Pentamer structures with var- 
ious proportions of HQ, IQ and SQ content have been constructed from 





Fig. 2. Schematic monomer structures used in this work and labelling of active 
sites: HQ- hydroquinone; IQ- 5,6-indolequinone; SQ- semiquinone. 
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geometrically optimized monomers and further analyzed for their spectro- 
scopic properties. All oligomer models have been built by using "Random 
Polymer Builder" module of the Cerius^/ MSI packet j23j . 

In assessing electronic spectral properties of modelled molecular struc- 
tures, we report here results from the spectroscopic ZINDO/s calculations 
developed by Zerner ^H] and co-workers. The method is based on re- 
parametrization of the semi-empirical INDO (Intermediate Neglect of Dif- 
ferential Overlap) approach that includes monoatomic differential overlap 
for one-centre integrals {i.e. integrals involving basis functions on the same 
atom) . These calculations yield absorption maxima which establish that the 
differences in mixing ratios of various monomeric units in an overall poly- 
meric structure of eumelanins can influence and alter absorption spectrum 
of the oligomer. 

2. Computational models 

In order to pinpoint the factors controlling the spectra of eumelanin 
polymers, the model subunits have been treated at several levels of struc- 
tural elaborations. As basic components of eumelanins 5,6-dihydroxyindole 
(HQ), 5,6-indolequinone (IQ) and semiquinone (SQ) (see Figure 2.) have 
been chosen. Since experimental data on the geometry of the molecules 
in Figure 2 are not available, the bare skeletons of indole molecules have 
been edited by a molecular editor in the Insightll/MSI module and fur- 
ther geometrically optimized by use of either ab initio Gaussian-94 or 
semiempirical consistent valence forcefield (CVFF) method. In both cases 
the energies for conformation were minimized until the convergence was 
achieved (the maximum derivative was less than 0.0001 kcal mor^l"^). 
The augmented CVFF program has been developed for material science ap- 
plications and is considered a standard forcefield to be used for organic 
molecules and polymers. As a default forcefield in Discover 95.0 molecular 
dynamics simulation package, it has been extensively used over the past 
years and is usually considered [201 12210^ well tested and characterized. It 
is primarily intended for studies of structures and binding energies, although 
is known also to predict satisfactorily frequencies in vibrational spectra as 
well as conformational energies. 

Structures optimized by CVFF method have been further used as input 
parameters for spectroscopic ZINDO |18| calculations. The method cir- 
cumvents cumbersome and time-consuming ab initio treatment of complex 
systems and is one of commonly used semiempirical molecular orbital (MO) 
methods. Over the last decade it has been succesfully applied in the study 
of bulk solids and defects in semiconductors . The INDO electronic cal- 
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culations are based on an all-valence-electron, self consistent field (SCF) 
molecular orbital procedure with the configuration interaction (CI). Within 
the general MO method, the Fock matrix can be represented as 

^ 1 

F^,u = V + E ^A,a[(A^J^|Aa) - -{^x\\ya)] (2.1) 

where h^j^y stand for the core integral 

V = / <j>^,{l)HcoreM'^)dTl, (2.2) 

(//i/|A(t), (/iA|z^cr) represent two-electron Coulomb integral and the exchange 
integral, respectively, 

{lJiv\\a) = <l)^il)M^)r^iM'^)M'^)dTidT2 (2.3) 

The electron density matrix P\u 

occ 

Pa. = 2^4,0,,, (2.4) 

i=l 

contains molecular orbital expansion coefficients c\i of the orbital formed 
as a linear combination of basis functions (f)\. The variation principle 
= applied to a closed-shell system leads then to the Roothaan-Hall 
equations: 

TV 

Y.{F^u-e^Sf,,)ci, = (2.5) 

for the set of orbital energies ej and MO coefficients Ci^. The elements of 
the matrix S^i, are the AO overlap integrals / 0* (l)0i,(l)(iT. 

Equation 1)2. 5|) is solved iteratively until a self-consistency is reached 
within the demanded accuracy. The final SCF solution yields desired MOs 
$j and their orbital energies e^. The ground state electron configuration is 
produced by filling the orbitals with all electrons in the order of increasing 
energy. 

The SCF calculations done with a set of M basis functions require the 
computation of matrix elements. In order to make a treatment of large 
molecules possible, one has to reduce the complexity by either replacing the 
effect of inner core electrons by effective (pseudo-) potentials or by applying 
the semi-empirical SCF method which neglects most of the matrix elements 
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(nvlXa) and parametrize the remaining ones. The term "intermediate ne- 
glect" points to the retention of one-center-two-electron exchange integrals 
{^v\lJLv). Within the Zerner's version of INDO (ZINDO) approximation jl8j . 
basis orbitals (pi are envisioned to be strongly orthogonal and they obey the 
relations 

{^,\B\\('a^) ^ 5ab5cd 1 1 <(l)</.?(l)rr2V?(2)</'f (2)dTidr2 
= (^^z^^lA^fj^) A = C 

= {fi^iy^\X^a^) A^C (2.6) 

where 6^ is the atomic orbital centered on atom A. In order to maintain ro- 
tational invariance, two-center integrals {A ^ C) are evaluated over atomic 
orbitals 4>^ that are s symmetric but have the same exponents and expansion 
coefficients as c;/)^. The one-center core integrals /i^^ are obtained from ion- 
ization potentials and the resonance integrals (the two-center h^i, integrals) 
are purely empirical parameters set to reproduce experimental spectra. The 
method as employed herein leads to the ground electronic states obtained as 
closed-shell molecular orbital wave functions in the restricted Hartree-Fock 
(RHF) framework. In the next step, low lying excited states were approx- 
imated by configuration interactions (CI) among configurations generated 
as single excitations from RHF ground state. The CI method included the 
highest 15 occupied and lowest 15 unoccupied molecular orbitals. ^ 

All polymers were arbitrarily chosen as pentamers built by use of the "Ran- 
dom Polymer Builder" module of the MSI packet |23j . Four different groups 
of polymers: homogenous polyHQ, polySQ, polylQ (c/. Figure 3) and het- 
eropolymers ( i.e. pentamers built up from basic units with random ratios 
of molecules HQ:IQ:SQ) were investigated. The polymer structures were 
formed by defining all possible pairs of "head" and "tail" positions out of 
5 labeled centers (see Figure 2) located on a monomer. Structural config- 
uration models of pentamers have been then developed along the polymer- 
ization direction in which "heads" were joined with "tails" of subsequent 
monomers. Resulting pentamer structures (c/. Figure 3) can be either lin- 
ear and flat, nonplanar linear or forked structures, respectively. Despite 
their various conformational structures, all three groups are characterized 
by a mean dimension (measured as a maximum observable distance between 
non-hydrogen atoms in a model pentamer form) of about 2nm and molec- 
ular weight about 720-740, i.e. with MW < 1000. If assuming that such 

^ Expanding the active basis of the CI space have not resuhed in significant differences 
in modelled absorption spectra, therefore all oligomeric structures have been treated 
within the same Cl-conformity. 
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units are basic chromophores in the eumelanin structure, the models predict 
that there are about of milhon of such pentamers in a melanin subunit of 
an average diameter 150 — 200nm [241 115j . 

Heteropolymeric structures have been created either with equal probabili- 
ties of choosing IQ, SQ or HQ molecule at each step of polymerization, or 
by using IQ and HQ monomers only, i.e with the HQ:IQ:SQ probability 
proportion taken as 3:1:0. The latter choice switches off the reduced form 
of IQ, the component SQ, that is mainly responsible for formation of free 
radicals produced during melanin synthesis P]. The planar and non-planar 
conformations were analyzed seeking for their effect on absorption spec- 
tra. Models of planar linear polymers (each monomer had no more than 2 
bonds), planar forked polymers and non-planar linear polymers (see Figure 
3.) have been investigated for their absorption spectra, calculated solely for 
neutral structures by using ZINDO-CI method. 

Dynamic interactions between the model basic "molecule" and the sur- 




Fig. 3. Examples of polymers used in calculations: A - planar and forked, B - 
planar and linear, C - nonplanar, linear. 
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rounding medium (composed of the set of randomly distribured similar 
oligomeric units) can be a source of fluctuations in the electronic gap energy 
between the ground and excited state. To account for these fluctuations and 
their effect on the system response to a radiation field, we have presented 
calculated spectra adopting a Brownian oscillator model (25j that relates 
variability in the electronic gap to displacements in some dynamical vari- 
able X representing the system under consideration. In the linear regime 
(weak coupling to the perturbing external radiation field), the fiuctuation- 
dissipation theorem relates the system's response function to the derivative 
of the correlation function C{t) =< x{t)x{t') > with brackets representing 
an equilibrium average. In the Brownian oscillator model the equilibrium 
states are defined by incorporating a model of the bath represented by a 
Gaussian stochastic force f{t) acting on a variable x: 

miXi{t) + miujfxi{t) + mi I drji^t - T)xj(r) = fi{t) + Fi{t) (2.7) 

where nii, Fi{t) stand for the "mass" of a mode Xi and an external driving 
force, respectively. The correlation of the random force is related to the 
dissipative term (friction), i.e. < fi{t)fk{t') >= 5ik2mikBT'yi{t — t'). The 
variation of the friction term with time (or frequency) reflects the timescale 
of the thermal motions of the bath. By assuming that the latter are very fast 
compared to the oscillator motion, we accept the ohmic dissipation model 
with 7j being a constant. In a strongly overdamped case (7^ >> 2uJi) and 
in a fast modulation limit [^Hj, the absorption lineshape calculated as 



with 



1 

(^absi^) = -Re dtex.p[i{uj - ujeg)t]ex-p[-g{t)] (2.8) 
TT Jo 

g{t) = f dt2 r dtiC{ti) (2.9) 

JQ Jo 

assumes a Lorentzian form 

r 

7r[(u; - UegY + r^J 

with the line width F proportional to the product of characteristic nuclear 
time scale (characteristic time scale of medium relaxation) and a coupling 
strength of the nuclear degrees of freedom to the electronic transition. With- 
out discussing the origin of the particular choice of the lineshape in more de- 
tails, we have thus arbitrarily adopted that the energy transitions calculated 
for the model system were Lorentzian enveloped and summed, weighted with 
an oscillator strength to simulate the overall absorption spectra. Differently 
chosen values of T reflect then the degree of dissipation of absorbed energy 
via coupling to the bath. 
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3. RESULTS 
3.1. Monomers 

Geometrical features of the ground state monomers used in this study 
have been compared with optimized coordinates from ab initio Gaussian 94 
calculations. In accordance with the previous analysis by Bolivar-Marinez 
et al. jl4j . various methods of optimization result in small variation of main 
bond lengths and bond and dihedral angles yielding structures that are basi- 
cally planar. This feature is further reflected in the major vr— vr* character of 
the HOMO-LUMO transition for all monomer molecules. To assess the reli- 
ability of our calculations based on CVFF method combined with ZINDO/s 
calculations in predicting spectral features of melanin's derivatives, we have 
calculated transition energies for the bare skeleton of HQ (with each omit- 
ted substituent replaced by a hydrogen atom), after optimization by either 
CVFF or PM3 method which has been used in a previous study ^3]. The 
first transition line calculated in ZINDO was 306. 5nm (296. 6nm) for the 
structure optimized with CVFF (PM3) program, respectively. This has 
been compared with the first transition line (300. 3nm) calculated for the 
" empirical" bare skeleton structure of tryptophan obtained from the Protein 
Data Bank. 



monomer 


Gaussian/ZINDO 


CVFF/ZINDO 




energy (osc. strength) 


energy (osc. strength) 




[nm] 


[nm] 


HQ 


307.2(0.05) 


316.2(0.06) 




228.0(0.82) 


232.6(0.87) 


IQ 


464.5(0.12) 


529.9(0.13) 




218.9(1.03) 


225.3(0.52) 


SQ 


481.0(0.02) 


856.4(0.03) 




217.4(0.94) 


238.2(1.46) 



Table 1: Calculated excitation energies [nm] for optimized structures of melanin's 
precursors. First row: first transition, second row: strongest transition. Oscillator 
strength are given in paranthesis. 

In all these cases, the first transition lines were predicted with oscillator 
strengths of about 0.02. By inspection, all structures displayed also similar 
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character of the calculated CI coefficients and were comparable with the 
previous theoretical and experimental 17 works related to melanin's 
monomers. 

Results of Gaussian/ZINDO and CVFF/ZINDO calculations on neutral 
HQ, IQ and SQ monomers are presented in Table 2. The first optically 
active electronic transitions (first row entry in Table 2) and the strongest 
active electronic transition (second row entry in Table 2) exhibit shifts to- 
wards longer wavelengths after CVFF structure optimization which allows 
structures with some degree of non-planarity. 

Reported transitions are dominated by a few configurations involving fron- 
tier orbitals. Inspection of the calculated CI coefficients reveals that about 
70% of the ffist excited state for HQ monomer is accounted for by HOMO^LU- 
MO excitation and about 15% of HOMO-l^LUMO+1 excitation. The 
strongest transition is composed of mixing among 2 highest occupied and 3 
lowest unoccupied orbitals. The pattern of the ffist excited state is essen- 
tially preserved for SQ and IQ monomers where >96% of excitation comes 
from HOMO— >LUMO transition with remaining excitations contributing 
with less than 1%. 



Figure 4 presents simulated optical absorption spectra of model poly- 
mers. The curves have been obtained from the optical transition spectrum 
enveloped by normalized Lorentzians weighted by the calculated oscillator 
strength; i.e energies (Aj) and oscillator strengths (sj) obtained from calcu- 
lations were converted to: 



with the width of a line set up arbitrarily to 10 (in units of A). The graphs 
refer to spectra of different mixtures of polymers: the first column displays 
spectra of the uniform mixture of polyHQ polymers {i.e. pentamers com- 
posed of HQ monomers, only), second and third columns show absorbance 
spectra of similar uniform mixtures of polylQ and polySQ, respectively. 
The fourth column presents spectra of heteropolymers composed with equal 
probabilities of choosing either one of three various monomers (i.e. with 
HQ:IQ:SQ=1:1:1) at each step of the "polymerization process" leading to a 
basic pentamer structure. The last column displays spectra of nonuniform 
polymers built by use of the proportion HQ:IQ:SQ = 3:1:0 at the basic 
level of modelling. Rows in Figure 4 relate to spectra obtained for different 
conformations of polymers. The first row represents mixtures of linear and 



3.2. Oligomers and Polymers 




(3.1) 
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planar polymers (c/. Figure 3), the second and third row display spectra of 
linear nonplanar and planar forked polymers, respectively. 

All spectra up to the last row in Figure 4 have been calculated as mix- 
tures of 15 pentamers, whereas the bottom row shows sums of simulated 
spectra for 45 polymers whose absorbance spectra are displayed in upper 
three rows of each column. As an example, the third column and second 
row of Figure 4 displays spectrum of 15 linear, nonplanar pentamers com- 
posed of SQ monomers, only. Accordingly, the inlet in 5th column and 3rd 
row presents absorption spectrum of 15 "forking" pentamers composed of 
HQ and IQ molecules taken in proportion 3:1 (with no contribution of SQ 
molecules). Sum of spectra for various conformational structures of that 
type is displayed in the bottom right corner of Figure 4. 
As can be inferred from Figure 4, spectra of polyHQ have dominant con- 
tributions within the range 200 — 400nm, whereas admixture of IQ or SQ 
components results in a red shift of the absorption lines. Conformation 
variations seem to result in varying intensity of absorption, it is ambiguous 
however, to conclude on the overall effect of conformations on spectra of 
model polymers. 

Figure 5 comprises model absorption spectra of "higher order" poly- 
mers: in the top left corner spectra of 45 polyHQ and 45 polylQ oligomers 
are displayed, second position in the first column represents spectra of 45 
polyHQ and 45 polySQ, whereas the bottom of first column contains sim- 
ulated spectra for mixture of 45 polylQ and 45 polySQ polymers. For 
comparison, spectra of mixture of three uniform oligomers (45 polyHQ + 
45 polylQ + 45 polySQ) have been displayed (second column and first row 
in Figure 5) along with 45 heteropolymers HQ:IQ:SQ=1:1:1 (second row) 
and 45 nonuniform oligomers HQ:IQ:SQ = 3:1:0 (second column, third row). 
Note that all "higher order" spectra of polymers are constructed as direct 
sums (mixtures) of spectra of pentamers and as such they do not take into 
account interactions among those units. 

Simulated optical absorption spectra suggest that the main contribution to 
the absorption in the range of 500 — 800nm may be attributed to the pres- 
ence of SQ molecules in oligomeric structures. Shape of simulated spectrum 
of mixtures of nonuniform polymers resembles experimental spectrum (Fig- 
ure 1) for synthetic DOPA melanin. The maximum intensity of modelled 
spectrum is, however, shifted to about 300nm, whereas the peak inten- 
sity in the experimental spectrum is observed already at 200nm. It can 
be expected that additional modulation of spectrum may arise from some 
environmental effects. For example, presence of localized charges may be 
a source of electrochromic shift in the spectrum with red or blue shifts in 
transition energies determined by the sign and location of charges |261 127j . 
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Fig. 4. Theoretical absorption spectra for different mixtures of polymers: first 
column - polyHQ, second column - polylQ, third column - polySQ, fourth and fifth 
columns - heteropolymers (HQ:IQ:SQ = 1:1:1 and HQ:IQ:SQ = 3:1:0, respectively). 
Rows refer to different conformational forms: first row - spectra of 15 linear planar 
polymers, second row - 15 linear non planar polymers, third row - 15 forked planar 
polymers. Last row represents sums of spectra in each column (45 polymers). 
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by higher content of charged forms of IQ monomers in model structures. In 
particular, it has been demonstrated in the previous theoretical study by 
Bolivar-Marinez et al. that negative ions of IQ, SQ and HQ molecules 
(-1) begin to absorb at about l.OeV. This effect is responsible, however, for 
the red shift in simulated optical absorption spectra of negatively charged 
monomers. 

Calculated optical spectra smoothen also and shift when assuming stronger 
dissipation of excitation energies in the system, resulting in a broader width 
of lines (see Figure 6). In the case of heterogeneous oligomers models 
(HQ:IQ:SQ = 3:1:0) fourfold higher line width parameter than originally 
chosen leads to a single humped spectrum with a broad, structureless tail. 
The effect, however, does not alter significantly absorbance in the 200nm 
region. 

Synthetic eumelanin whose properties have been modelled in this study 
lacks a protein coat and therefore its optical spectrum does not match the 




200 
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Fig. 5. Simulated spectra for "higher order" mixtures of polymers: left column, top 
- 45 polyHQ + 45 polylQ, below - 45 polyHQ + 45 polySQ, bottom - 45 polylQ 
+ 45 polySQ. Right column - 45 polyHQ + 45 polylQ + 45 polySQ, below - 45 
heteropolymers (HQ:IQ:SQ = 1:1:1), bottom - 45 heteropolymers (HQ:IQ:SQ = 
3:1:0). See text for details. 
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natural pigment, especially in the region about 250 — 270nm, where existing 
evidence jll] of spectral features has been attributed to absorption of light 
by protein environment. 

4. CONCLUSIONS 

The absorption spectra for model of eumelanins were calculated as- 
suming melanin to be composed of mixture of various unit-polymers. All 
molecules were parametrized as random pentamers of HQ, SQ and/or IQ. 
The idea of "mixed nature" of the pigment, composed of heterogeneous 
oligomers in a random, not necessarily well organized assemblies is compat- 
ible with experimental data. Simulated optical spectra manifest the effect 
of aggregation in oligomer conformation resulting in broadening of the ab- 
sorbance arm in the region of 500 — 800nm. Peak intensity of spectra can 
be correlated with the presence of neutral HQ monomers in model struc- 
tures. At higher level of aggregation (135 pentamers and more), further 
addition of a polymeric unit (a pentamer) to the existing structure does not 
alter the character of the spectrum. Also, conformational deformations of 




200 500 800 1100 

Wavelength [nm] 

Fig. 6. Effect of assumed line width on spectra of 45 heteropolymers (built from HQ 
and IQ monomers only). Starting from the top, the values of Lorentzian's widths 
have been set up to 10, 20 and 40 (in units of the wavelength), respectively. 
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constituents (oligomers) have only a minor effect on the structure of mod- 
elled optical spectra. Our future studies aim at understanding stabilisation 
of radicals in a system of stacked oligomers and interpretation of factors 
responsible for radical-trapping processes in eumelanins. 
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